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LONG-TERM  GOALS 

The  long-term  goals  of  this  project  are  to  obtain  an  understanding  of  the  mechanisms  by  which  tidal 
energy  is  used  to  vertically  mix  the  ocean  against  the  action  of  gravity.  Ultimately  better 
parameterizations  of  the  mixing  caused  by  tides  will  result,  allowing  better  prediction  of  coastal 
dynamics,  biogeochemistry  and  sediment  transport  and  the  oceanic  general  circulation. 

OBJECTIVES 

The  process  of  mixing  by  tides  interacting  with  topography  involves  several  stages.  First  some 
fraction  of  the  energy  contained  in  the  barotropic  tide  must  be  converted  into  baroclinic  energy, 
through  the  generation  of  internal  tides  and  turbulent  boundary  layers.  Secondly,  the  energy  in  the 
internal  tides  must  be  transmitted  into  smaller  vertical  wavelengths,  thereby  increasing  the  vertical 
shear  of  the  motion.  When  vertical  shears  are  sufficiently  strong,  instability  may  result,  leading  to 
overturning  and  mixing.  Finally,  the  mixed  fluid  is  transported  away  from  the  mixing  region 
modifying  the  ocean  stratification.  The  net  effect  of  the  tides  on  the  ocean  stratification  depends  on  the 
efficiency  of  all  three  processes. 

Our  objectives  are  to  understand  (a)  the  generation  of  internal  tides  by  the  interaction  between 
barotropic  tides  and  topography  including  finite-amplitude  three-dimensional  variations  in  topography, 
finite-amplitude  barotropic  tidal  forcing,  non-hydrostatic  effects  and  the  boundary  layer  processes;  (b) 
the  mixing  generated  by  internal  tides  reflecting  from  a  sloping  boundary  in  the  presence  of  both  two- 
and  three-dimensional  variations  in  slope,  and  finite  rotation;  and  (c)  the  mechanisms  of  lateral  and 
isopycnal  transport  of  mixed  fluid  away  from  the  boundary  induced  by  the  secondary  circulations 
generated  through  spatial  variations  in  mixing.  Earlier  studies  have  ignored  three-dimensional  large 
amplitude  variations  in  topography  and  non-hydrostatic  effects  (which  are  important  for  small-aspect 
ratio  motion). 

APPROACH 

We  use  high-resolution  numerical  simulations  to  explicitly  resolve  the  turbulent  mixing  processes.  For 
such  simulations  we  require  a  numerical  model  which  can  (a)  capture  the  non-hydro  static  physics  of 
overturning  and  mixing  processes  and  (b)  include  arbitrary  three-dimensional  variations  in  topography. 
The  Marshall  et  al.  (1997a,b)  code  (known  as  the  MIT  ocean  model),  which  is  non-hydrostatic,  and 
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includes  topography  through  a  finite -volume  formulation,  is  such  a  model.  We  are  carrying  out  three 
different  groups  of  simulations:  (a)  We  impose  topography,  baro tropic  tides  and  subinertial  flows 
suggested  by  recent  observations  (e.g.  the  region  of  the  TWIST  field  program,  Polzin  et  al,  1998  and 
Monterey  Canyon  region,  Petruncio  et  ah,  1998),  and  investigate  the  internal  tide  generated  by  the 
flow-topography  interactions,  comparing  these  results  with  earlier  models  which  assume  small- 
amplitude  (e.g.  Bell,  1975)  or  two-dimensional  finite  amplitude  (Baines  1982)  topography,  (b)  We 
impose  internal  tide  forcing  and  investigate  the  interaction  with  a  two-dimensionally  varying  slope, 
focusing  on  the  influence  of  finite  rotation  on  the  overturning  and  mixing,  and  the  effect  of  slope 
variations  in  localizing  mixing,  comparing  with  earlier  laboratory  and  numerical  studies  with  uniform 
slope  and  no  rotation  (Ivey  and  Nokes,  1989;  Slinn  and  Riley,  1996).  (c)  We  examine  reflection  of 
internal  tides  from  three-dimensionally  varying  topography. 

WORK  COMPLETED 

In  the  past  year  we  have  focused  on  analysis  of  simulations  of  both  internal  tide  generation  and 
reflection  in  the  TWIST  region,  working  closely  with  Kurt  Polzin  to  compare  results  with 
observations,  and  simulations  of  internal  tide  reflection  and  breaking  on  two-dimensional  variable 
slopes.  We  have  carried  out  the  following  series  of  simulations. 

Data  from  the  LIWI  funded  TWIST  project  (Kurt  Polzin,  PI)  have  been  used  to  initialize  the 
stratification  and  topography  for  a  region  of  the  East  Coast.  The  model  has  been  forced  with  either  the 
cross-slope  barotropic  tidal  signal  (imposed  at  the  off-shore  open  boundary)  or  with  the  along-slope 
barotropic  tidal  signal  (imposed  through  a  body  forcing  term).  The  generation  of  internal  tides  at  the 
shelf-break,  and  over  the  corrugations  on  the  slope  is  simulated,  and  the  scattering  of  the  shelf-break 
waves  from  the  corrugations  examined. 

Several  two-dimensional  simulations  of  internal  tide  breaking  on  a  continental  slope  of  variable 
topography  have  been  completed. 

In  addition  to  these  simulations,  theoretical  work  has  been  completed  to  help  interpret  the  results, 
including  (a)  an  extension  of  Thorpe  (2001)  to  include  rotation,  necessary  to  describe  the  scattering  of 
the  shelf-break  internal  tide  from  corrugations  and  (b)  development  of  a  wave  Froude  number  criterion 
for  the  range  of  topographic  slopes  over  which  internal  wave  breaking  will  occur. 

Three  manuscripts  have  been  completed  (Legg  and  Polzin,  2002;  Legg  2002;  Legg  and  Adcroft,  2002) 
and  submitted  for  publication. 

Under  my  supervision,  a  summer  student  has  carried  out  a  series  of  two-dimensional  calculations  to 
investigate  the  role  of  nonlinearity  (large  barotropic  flow  speed)  in  modifying  the  barotropic  to 
baroclinic  energy  conversion,  and  the  ratio  of  mixing  to  waves  for  the  fate  of  that  energy,  for  both 
small  amplitude  and  large  amplitude  topography. 

RESULTS 

In  simulations  with  TWIST  topography  and  stratification  an  internal  tide  is  generated  at  the  shelf-break 
by  the  cross-slope  barotropic  tide.  This  internal  tide  reflects  from  the  slope  in  a  region  of  corrugations; 
in  addition  the  reflected  wave  scattered  waves  are  generated  with  an  along-slope  wavenumber  given  by 
the  wavenumber  of  the  corrugations.  These  scattered  waves  are  not  present  if  planetary  rotation  is 
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ignored,  and  result  from  the  flow  of  the  along-slope  velocity  component  of  the  primary  internal  tide 
(zero  in  the  absence  of  rotation)  over  the  corrugations.  The  scattered  waves  can  be  predicted  by 
extending  the  theory  of  Thorpe  (2001)  to  include  rotation,  and  there  is  qualitative  agreement  between 
the  forms  of  the  predicted  velocity  profiles  and  the  simulated  velocity  profiles,  in  addition  to  the 
observed  velocity  profiles.  Hence  wave  scattering  at  corrugated  topography  is  a  mechanism  for 
scattering  low  mode  internal  wave  energy  into  higher  vertical  modes,  generating  greater  shear  and 
greater  potential  for  mixing. 


(a)  2D  wiLh  rolALion 


(h)  3D  wiLh  roLAlion 


(c)  2D  no  rolALion 


(<d)  3D  no  roLALion 


Figure  1:  Vertical  profiles  of  cross-slope  baroclinic  velocity  at  a  location  on  the  continental  slope 
for  several  tidal  cycles,  shown  for  several  different  simulations:  (a)  2D  calculation  with  rotation,  (b) 
3D  calculation  with  rotation,  including  topographic  corrugations  of  wavelength  3.3km,  orientated 
up  and  down  the  slope,  (c)  2D  calculation  without  rotation  (d)  3D  calculation  without  rotation, 
including  corrugations.  When  both  rotation  and  corrugations  are  included,  much  more  small-scale 
structure  is  seen  in  the  velocity  field,  leading  to  greater  shear  especially  at  depth. 


Mixing  generated  by  internal  tide  reflection  at  concave  and  convex  slopes 

Our  numerical  calculations  have  examined  the  reflection  of  an  internal  wave  from  continental  slope 
topographies  which  are  planar,  concave  or  convex.  Whereas  earlier  analytical  studies  predict  that 
concave  slopes  near  the  critical  angle  will  lead  to  less  mixing  than  planar  or  convex  slopes,  we  find 
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strong  mixing  in  all  cases,  in  the  form  of  bores  propagating  up  the  slope,  as  well  as  shear  instability 
behind.  In  all  cases  a  region  of  mixed  fluid  results  bounded  below  by  the  topography,  and  above  by 
the  wave  characteristic.  We  conclude  that  the  principal  difference  between  our  calculations  and  the 
earlier  studies  is  the  nonlinearity  of  the  flow.  The  bores  result  when  the  Froude  number  (flow  speed  / 
wave  group  velocity)  of  the  reflected  wave  exceeds  unity  -  the  reflection  process  tends  to  increase  the 
velocity  amplitude  and  decrease  the  group  velocity  of  the  wave.  We  derive  a  range  of  slope  angles 
around  the  critical  angle  at  which  internal  tide  breaking  can  occur,  and  our  predictions  agree  well  with 
our  numerical  results. 


Snapshots  of  IJ  and  b  at  t=5  39  T 


(a)  f*tarxfzr 
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(b)  (Jorxcmix 


(c)  Corjuxx 


Figure  2:  Snapshots  of  cross-shelf  velocity  (color)  and  isopycnals,  from  3  different  simulations  of  a 
mode  1  internal  wave  reflecting  from  a  sloping  boundary.  In  all  three  cases  the  average  slope  is  at 
the  critical  angle  for  the  incident  internal  wave,  but  the  slope  shape  is  (a)  planar,  (b)  concave,  (c) 
and  convex  about  the  mid-depth  critical  point.  Whereas  analytical  theory  predicts  enhanced  mixing 
due  to  critical  angle  reflection  only  for  planar  and  convex  slopes,  all  three  simulations  show 
turbulent  mixing,  caused  by  an  internal  bore  propagating  up  the  slope.  The  bore  is  characterized  by 
sharp  horizontal  density  gradients,  and  local  overturning. 


IMPACT/APPLICATIONS 

Our  results  are  helping  in  the  interpretation  of  observations  of  tidally  forced  flows  on  the  continental 
slope  observed  by  LIWI  investigators  and  others.  Our  results  are  demonstrating  mechanisms  for 
exciting  small  vertical  scale  baroclinic  response  to  barotropic  flow  over  topography.  These  small 
vertical  scales  (and  hence  large  shears)  are  essential  if  the  tidal  forcing  is  to  generate  mixing,  as  widely 
speculated  in  recent  years.  Ultimately  these  results  may  help  to  improve  parameterizations  of  tidal 
mixing  over  topography  for  inclusion  in  global  and  coastal  circulation  models. 
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TRANSITIONS 


We  are  actively  communicating  results  with  K.  Polzin,  E.  Kunze,  J.  Nash,  R.  Street,  and  R-C  Lien,  to 
help  interpret  observations  made  during  LIWI  and  compare  with  other  numerical  simulations  of  these 
phenomena. 

I  am  currently  discussing  ways  of  implementing  the  results  of  the  internal  tide  breaking  study  in  a 
parameterization  of  boundary  mixing  with  researchers  at  NOAA  GFDL  and  members  of  the  MIT 
ocean  modeling  development  team. 

RELATED  PROJECTS 

This  project  examines  processes  closely  related  to  observations  included  in  LIWI  (Polzin,  Toole  and 
Schmitt  and  Paduan,  Rosenfeld,  Kunze  and  Gregg). 
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